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Executive  Summary 

We  are  developing  a  comprehensive  computational  and  experimental  program  to  investigate  the 
characteristics  and  dynamic  behavior  of  nanosize  aluminum  in  novel  energetic  materials.  Two 
classes  of  advanced  propellants  will  be  considered.  The  first  will  be  composite  propellants,  consisting 
of  solids  including  AP,  HMX,  RDX,  or  nanoaluminum  with  a  suitable  binder  such  as  HTPB.  The 
second  will  be  a  nanoscale  aluminum  and  ice  propellant.  For  the  computational  portion,  we  will  (i) 
use  our  in-house  packing  code  to  generate  morphologies  of  interest  for  heterogeneous  propellants; 
(ii)  modify  our  CFD  code  to  include  appropriate  models  for  chemistry,  radiation,  and  nano-sized 
aluminum;  and  (iii)  simulate  transient  behavior  of  the  propellants  at  rocket  motor  conditions.  For 
the  experimental  portion,  we  will  (i)  use  the  mixing  facility  at  Purdue  to  prepare  appropriate 
propellant  samples  containing  nano-sized  aluminum,  thus  allowing  for  a  full  characterization  of 
the  propellant  (composition;  particle  size  distribution)  necessary  for  a  coordinated  modeling  and 
experimental  program;  and  (ii)  carry  out  both  steady  and  unsteady  experiments  to  fully  characterize 
the  propellants.  The  experiments  will  be  carried  out  in  a  window  pressure  vessel  rated  to  6000  psi 
and  in  a  small  end-burning  rocket  motor.  The  experiments  will  also  be  used  for  model  calibration 
and  for  validation. 

This  final  report  describes  Phase  1  work  on  the  computational  modeling,  validation,  and  pro¬ 
pellant  formulation.  Initial  results  show  excellent  agreement  with  prior  experimental  data,  as  well 
as  the  path  forward  for  the  next  steps  of  the  work. 

1  Introduction 

This  final  report  documents  the  status  of  initial  work  performed  on  STTR  contract  FA9550-09-C- 
0012,  Dynamic  Behavior  of  Nano- Sized  Particles  in  Novel  Energetic  Materials  for  Space  Propulsion. 
The  outline  of  the  report  is  as  follows.  We  begin  by  restating  the  Phase  I  Technical  Objectives  and 
the  Phase  I  Work  Plan.  In  Section  2  we  present  a  literature  review  of  relevant  work.  In  Section  3  we 
present,  in  its  entirety,  the  mathematical  model  for  nano-aluminum-based  propellants.  Although  a 
portion  of  the  model  appears  in  [1],  we  felt  it  best  to  present  a  clear  picture  for  future  reference. 
Results  are  presented  in  Section  4.  As  part  of  Phase  I  funding,  we  performed  a  validation  study 
in  Section  4.1  by  comparing  burn  rates  from  the  model  and  simulations  with  those  obtained  from 
the  experimental  data  of  2  .  The  effort  involved  generating  a  model  for  the  morphology,  modifying 
the  radiation/combustion  code,  performing  a  grid  resolution  study,  and  then  simulating  the  burn 
through  of  the  model  propellant  for  various  pressures.  We  note  that  excellent  agreement  is  obtained. 
In  Section  4.2  we  present  preliminary  results  for  a  propellant  manufactured  at  Purdue  as  part  of 
Phase  I  funding.  A  more  detailed  presentation  of  the  Purdue  experimental  results  are  given  in 
Section  5. 

1.1  Phase  I  Technical  Objectives 

Our  research  and  development  program  investigated  the  applicability  of  our  simulation  framework 
for  novel  energetic  higli-density  propellants,  established  the  feasibility  of  an  experimental  validation 
program  for  these  propellants,  and  provided  a  foundation  for  commercializing  the  simulation  soft¬ 
ware.  The  simulations  were  performed  at  IllinoisRocstar  LLC  and  the  experiments  were  performed 
at  Purdue.  To  this  end  our  Phase  I  objectives  were  as  follows: 
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1.  Develop  a  matrix  of  propellant  formulations  by  varying  AP  particle  size,  solids  content  (AP, 
HMX,  RDX,  nano-aluminum),  pressure,  initial  propellant  temperature; 

2.  Prepare  samples  with  the  mixing  facility  at  Purdue; 

3.  Conduct  preliminary  burning  rate  experiments  at  pressures  up  to  0000  psi,  measuring  burning 
rate  and  using  high-speed  imaging  and  high-speed  spectroscopy; 

4.  Conduct  preliminary  dynamic  experiments  of  ignition  under  realistic  rocket  conditions  in  a 
windowed  pressure  vessel  by  a  CO  2  laser; 

5.  Modify  the  current  simulation  code  to  include  models  for  nitramines,  nano-aluminum,  radi¬ 
ation,  and  combustion.  As  discussed  above,  there  is  already  in  existence  a  substantial  foun¬ 
dation  for  the  work,  but  certain  modifications  will  have  to  be  made  for  the  novel  energetic 
propellants  proposed  here; 

A  number  of  these  objectives  were  met  in  Phase  1. 

1.2  Phase  I  Work  Plan 

Under  Phase  I  STTR  funding  we  firmly  established  the  feasibility  of  using  our  simulation  codes  for 
the  new  class  of  propellant  formulations  proposed  here.  We  also  developed  the  experimental  valida¬ 
tion  program,  and  planned  the  commercialization  process  for  the  engineering  analysis  product (s). 
In  the  remainder  of  this  section,  we  will  expand  the  description  of  the  technical  issues  introduced 
in  the  previous  section. 

1.  The  computational  approach  employed  our  in-house  packing  code  to  create  three-dimensional 
random  packs  of  spheres  that  can  represent  the  oxidizer  and  nitramines.  The  nano-sized 
aluminum  is  treated  by  homogenizing  it  with  the  binder  to  create  a  blend.  The  thermo- 
mechanical  behavior  (density,  thermal  conductivity,  regression  rate)  is  modeled  using  our 
homogenization  strategies.  The  energetic  consequences  of  the  addition  of  nano-sized  alu¬ 
minum  in  the  gas-phase  combustion  field  is  modeled  by  a  radiation  transport  equation  and 
by  aluminum  oxidation  heat  release. 

2.  The  experimental  approach,  carried  out  at  Purdue,  focused  on  both  steady  and  unsteady 
behavior  of  novel  energetic  propellants.  After  a  complete  matrix  has  been  properly  defined 
(AP,  HMX,  RDX,  nano-aluminum,  HTPB:  with  input  from  the  sponsor  or  AFRL  personnel), 
some  samples  were  created  by  the  mixing  facility  at  Purdue  as  part  of  an  overall  feasibility 
study.  Each  sample  is  characterized  with  respect  to  composition  and  particle  size  distribu¬ 
tion.  This  allowed  for  a  coordinated  effort  between  the  modeling  and  the  experiments.  The 
samples  is  then  placed  in  the  windowed  pressure  vessel,  where  both  the  pressure  and  the 
initial  temperature  of  the  sample  is  varied.  We  note  that  most  previously  reported  work  on 
propellants  focus  mainly  on  pressure;  here,  we  varied  the  initial  temperature  as  part  of  the 
study  matrix.  The  burning  rate,  high-speed  imaging,  and  high-speed  spectroscopy  is  obtained 
for  the  prepared  samples. 
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2  Literature  Review 


Many  factors  are  important  when  selecting  propellant  formulation  in  a  solid  motor  system.  Of 
special  importance  is  the  notion  of  internal  ballistics,  a  category  that  encompasses  propellant  char¬ 
acteristics  such  as  burn  rate,  physical  properties,  and  burning  surface  phenomena.  Furthermore, 
oxidizer,  fuel,  and  binder  selection  will  also  affect  other  important  design  decisions  such  as  specific 
impulse,  ignitability,  and  aging.  In  the  strive  to  maximize  Isp  and  increase  propellant  burn  rates, 
designers  have  typically  resorted  to  energetic  catalysts,  increased  oxidizer  loading,  or  reducing  the 
oxidizer  particle  size  [3].  In  ammonium  perchlorate,  aluminum,  and  hydroxyl- terminated  polybuta¬ 
diene  (AP/A1/HTPB)  composite  propellants,  the  primary  topic  of  interest  for  this  project,  further 
optimization  has  also  resulted  in  the  use  of  bimodal  coarse  and  fine  AP  distributions  to  modify 
flame  kinematics  and  improve  surface  burning  where  the  introduction  of  fine  oxidizers  generally  re¬ 
sult  in  an  increase  in  burning  rate  [5].  Conversely,  studies  of  propellant  with  decreasing  micrometric 
aluminum  sizes  show  only  a  marginal  increase  [14].  When  transitioning  to  nano-scale  aluminum 
particles,  though,  the  drastically  increased  reactive  surface  area  per  unit  mass  results  in  completely 
different  burning  characteristics  leading  to  reduced  agglomeration,  higher  combustion  temperature, 
and  lower  ignition  delay  [10],  warranting  further  investigation.  To  accomplish  a  characterization  of 
the  dynamic  behavior  of  AP/nAl  composites,  transient  effects  will  be  studied  in  Crawford  bomb 
strand  burners  and  a  sub-scale  rocket  motor. 

Before  pursuing  such  endeavors  it  is  first  necessary  to  formulate  an  appropriate  propellant  test 
matrix  based  on  existing  work  [2-5].  Doing  so  will  allow  verification  of  Purdue’s  mixing  procedures 
and  repeatability  tests,  which  will  then  allow  robust  calibration  for  the  CFD  code  being  developed 
at  IllinoisRocstar,  LLC.  For  completion,  a  cursory  overview  of  considerations  in  AP/A1  formulation 
selection  is  presented  in  sections  2.1  and  2.2. 

2.1  Ammonium  Perchlorate 

As  stated  earlier,  the  burn  rate  of  a  solid  propellant  can  be  increased  by  increasing  the  amount  of 
loaded  oxidizer,  reducing  the  oxidizer  particle  size,  or  both.  However,  this  can  only  go  so  far  before 
manufacturing  considerations  put  a  constraint  on  maximum  amount  of  allowable  solids  (typically 
anything  higher  than  90%  solids  cannot  be  mixed).  To  understand  the  role  of  AP  size  in  propellants 
it  is  necessary  to  introduce  the  role  of  the  leading  edge  flame  (LEF)  discussed  by  [9].  As  an  example, 
[5]  compared  the  burn  rate  effects  of  AP  and  aluminum  formulation  at  pressures  up  to  5000  psi 
where  they  concluded  that  direct  replacement  of  pAl  with  nAl  shows  virtually  no  change  when 
using  monomodal  AP.  By  contrast  bimodal  AP  exhibited  markedly  different  behavior  over  the 
entire  pressure  range. 

Related  to  this  is  the  fi A1  agglomeration  sizes  as  studied  by  [4j.  Here  the  authors  discuss  how 
large,  monomodal  AP  creates  oxidizer-rich  regions  which  in  turn  results  in  unfavorable  ignition 
conditions.  Consequently  the  aluminum  liquefies  and  accumulates  before  reacting  resulting  in 
increased  slag  formation  and  unburned  fuel.  The  specifics  however,  are  dependent  on  the  ratio  of 
fine  to  coarse  AP  as  well  as  the  selected  sizes,  but  both  [5]  and  [4]  illustrate  the  need  for  bimodal  AP 
to  advantageously  utilize  nano-aluminum  as  an  energetic  replacement  for  micrometric-aluminum. 
In  general,  though,  burn  rate  is  most  enhanced  when  the  ratio  of  fAP  to  cAP  is  approximately  0.25 

[ii]- 
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2.2  Nano- Aluminum 

Many  theories  have  been  presented  to  explain  the  marked  increase  in  burn  rate  when  using  nAl 
including  the  following: 

1.  Shortened  ignition  delay  since  less  time  is  needed  to  heat  a  nanoparticle  to  ignition  temper¬ 
ature  2]; 

2.  Burning  at  the  surface  of  the  propellant  results  in  additional  conductive  and  radiative  heat 
interaction  [12]; 

3.  Thermal  expansion  of  nano-aluininum  is  enough  to  crack  the  noil-reactive  aluminum  oxide 
coating,  exposing  unreacted  aluminum  [10: ; 

4.  Increased  specific  surface  area  results  in  more  readily  reactive  surface  atoms  [10]. 

No  matter  the  reason,  the  increased  reactivity  of  nano-aluminum  compared  to  micrometric 
aluminum  requires  additional  attention  to  detail  regarding  storage,  mixing,  and  ingredient  compat¬ 
ibility.  Despite  the  many  types  and  origins  of  nAl  available,  only  three  will  be  considered  due  to 
their  immediate  availability  at  Purdue,  as  summarized  in  Table  1. 


Table  1:  Nano- Aluminum  Inventory  at  Purdue 


Manufacturer 

Origin 

Designation 

Nominal  Size  [nm 

Comments 

Argouide 

Russia 

ALEX 

50 

Uncoated 

Argon  ide 

Russia 

L-ALEX 

50 

Palmitic  Acid  Coating 

NovaCentrix 

USA 

N80 

80 

Uncoated 

This  recently  acquired  stock  has  not  yet  undergone  SEM  imaging  for  size  distribution,  BET 
analysis  to  determine  the  specific  surface  area,  or  hydrochloric  acid  decomposition  to  determine 
the  percent  active  aluminum,  but  historical  data  in  Table  2  indicates  that  the  palmitic  acid  coating 
greatly  increases  the  percent  active  aluminum,  the  ratio  between  aluminum  and  aluminum  oxide, 
AI2O3,  while  the  ALEX  and  N80  were  noticeably  lower.  As  emphasized  by  [7]  the  storage  of  un¬ 
coated  nAl  is  extremely  important  as  they  were  able  to  demonstrate  a  reduction  of  active  aluminum 
in  ALEX  from  80%  to  30%  in  seven  days  at  accelerated  conditions  of  40° C  and  75%  relative  hu¬ 
midity,  and  although  Purdue  raw  stock  is  stored  in  an  inert  argon  environment,  cast  propellant  is 
stored  at  ambient  conditions  and  may  be  subject  to  time-dependent  aging  and  moisture  absorption. 
In  fact,  a  study  using  100-150nm  ALEX  and  100-150nm  C-ALEX  (viton-A  coated)  showed  that 
C-ALEX  increased  burn  rate  126%  over  pure  HTPB  solid  fuel,  whereas  the  ALEX  only  yielded  a 
61%  increase.  Consequently,  unnecessary  exposure  of  either  raw  or  cast  nAl  to  the  atmosphere  is 
necessary  to  ensure  repeatable  mixes  and  results. 

In  addition,  the  high  specific  surface  area,  s,  compared  to  [2]  whose  nAl  typically  had  s  < 
15  m2/g,  poses  potential  mixing  problems  since  the  viscosity  of  the  uncured  propellant  increases 
drastically  wnth  s.  Also,  the  high  specific  surface  area  also  increases  the  nAl’s  likelihood  of  spon¬ 
taneous  reaction  with  either  ambient  moisture  or  oxidizing  agents  during  casting.  Both  issues  are 
documented  in  section  4. 
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Table  2:  Historical  nAl  Characteristics  (October  2007) 


Designation 

Average  Particle  Size 

nm 

Specific  Surface  Area  [m2/g] 

Percent  Active  Aluminum 

ALEX 

n/a 

n/a 

74 

L-ALEX 

n/a 

n/a 

90 

N80 

80 

25.7 

77.7 

3  Mathematical  Model 


3.1  Radiative  Transfer  Equations 

The  cloud  formed  by  the  fine/ultrafine  particles  is  modeled,  from  the  radiative  standpoint,  as  a 
gray  medium  with  time-space  varying  emission  and  absorption  characteristics.  In  this  section  the 
radiative  heat  transfer  equation  for  the  cloud  will  be  derived.  The  derivation  assumes  that  the 
propellant  surface  lies  in  the  x  —  y  plane  and  burns  in  the  nominal  ^-direction. 

The  radiative  transfer  equation  (RTE)  for  a  gray  medium  is  given  by 


S-VI  +  (Ka  +  Ka)I  =  KemIh  +  ^  [  I*dn\ 

4tt  j47r 

v  =  ( ,  — 

V  Ox '  dy '  Oz 

r  r‘2n  rTT 

dll' =  /  sin  W#', 

Jin  Jd>'=0  Jfl'= 0 


(1) 

(2) 

(3) 


where  I  =  /(r,  .s )  is  the  intensity  that  depends  on  position  and  direction,  <1>  =  <h(s,  s')  the  scattering 
phase  function,  r  =  (x,y,z)  the  position  vector,  and  s  =  (sx,sy,sz)  the  direction  cosines 


,sx  =  sin  0  cos  <fi,  sy  —  sin  0  sin  0,  sz  =  cos  0 , 


(4) 


with  0  <  0  <  7r  the  polar  angle  with  respect  to  the  z-axis,  0  <  (f>  <  2n  the  azimuthal  angle  with 
respect  in  the  x  —  y  plane,  and  the  blaekbodv  intensity  h  is  given  by 

Ih  =  aTA/n ,  (5) 

with  a  =  5.67  x  10-12  W/cm2K4,  the  Stefan-Boltzmann  constant.  The  preceding  parameters  are 
the  scattering  coefficient  Ks,  the  absorption  coefficient  Ka,  and  the  total  emission  coefficient  Kem. 
The  incident  radiation  G  and  the  radiative  heat  flux  vector  q  are  defined  as 

G(r)  =  [  I{r,s')dtt\  (6) 

q(r)  =  f  I{r,s’)s'd{Y.  (7) 

/At: 

Our  model  makes  a  number  of  assumptions: 

1.  The  particulate-mixture  cloud  can  be  modeled  as  a  gray  medium;  the  radiative  properties 
(absorption  and  scattering)  are  independent  of  the  spectral  wave  number. 
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2.  The  propellant  is  opaque  (i.e.,  the  transmissivity  r  is  zero). 

3.  We  ignore  in  depth  absorption  at  the  surface. 

4.  The  propellant  surface  is  a  ditfuse-gray  surface  and,  as  such,  it  emits  and  reflects  diffusely 
(i.e.,  the  exiting  intensity  is  independent  of  direction);  then 

I(rw,  s)  =  Jw{rw)/n,  n  •  s  >  0,  (8) 

where 

Jw  =  a rlb  4-  pH  (9) 

is  the  surface  radiosity  and  the  subscript  w  denotes  a  value  on  the  propellant  surface.  Here, 
H  is  the  hemispherical  irradiation  (i.e.,  incoming  radiative  heat  ttux) 

H(rw)=  (  J(ru„s)|n  •  s'\dSY,  (10) 

J  n-8r<  0 

where  n  is  the  local  outward  normal,  and  n  s'  =  cos#'  is  the  cosine  of  the  angle  between 
any  incoming  direction  s'  and  the  surface  normal. 

5.  The  intensity  is  periodic  in  the  x  and  y  directions  (the  propellant  burns  downward  in  the 
nominal  — z  direction). 

6.  The  phase  scattering  function  is  taken  to  be  linear-anisotropic 


4>(s,  s)  =  1  -f  A\s  •  s'.  (11) 

Setting  A\  —  0  implies  isentropic  scattering.  However,  Brewster  and  Parry  [15]  note  that 
smoke  is  anisotropic  and  so  we  allow  for  this  possibility. 

Most,  if  not  all,  of  the  assumptions  could  be  relaxed  if  more  information  about  the  radiation 
field  is  known.  In  any  event,  we  wish  to  begin  the  model  within  the  simplest  possible  framework. 

To  solve  the  radiative  transfer  equation  (1)  we  use  the  P\  approximation  of  the  spherical  har¬ 
monics  method,  which  is  equivalent  to  the  first  moment  method.  We  begin  by  writing  for  the 
intensity 

/(r,  a)  =  a(r)  -f  b(r)  ■  s.  (12) 

Substituting  (12)  into  (6)  and  (7)  yields  the  relations 

47T 

G  =  47T  a,  q  =  —b.  (13) 

which  in  turn  yields  for  the  intensity 

=  ^\G{r) +  3q{r)  ■  s\.  (14) 

In  deriving  these  relations  we  used  the  following  identities 

j  dCl1  =  4ir,  f  s’diV  =  0,  /  s’s'diV  =  —1,  (15) 

J  4n  J  J  4n  3 
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where  T  is  the  unit  tensor.  The  integral  on  the  right  hand  side  of  the  RTE  (1)  then  becomes 


J(r,  *')*(* 


G  4-  A\q  •  s . 


(16) 


Substituting  the  intensity  (14)  into  (1),  multiplying  by  (l,s),  respectively,  and  integrating  each 
equation  over  all  solid  angles,  yields  the  following  two  relations 


V  •  q  =  AvKemh  -  KaG  =:  /?,  (17) 

VG  =  -(/Vs  +  Ka)[ 3  -  iM*  =:  P</,  (18) 

where  uj  =  KS/(KS  4  Ka)  is  the  albedo.  These  last  two  relations  can  be  combined  to  yield  a  single 
equation  for  G  in  conservative  form: 


V2C 


VP  - VC 
p 


=  PR. 


(19) 


This  equation,  subject  to  appropriate  boundary  conditions  that  will  be  defined  below,  is  to  be 
solved  numerically  in  the  gas-phase  together  with  the  zero  Mach  number  reactive  equations.  Once 
G  is  known,  the  radiative  heat  flux  is  given  by 


and  in  the  energy  equation  we  have 


(20) 


V  q  —  47rA(:m/fr  KaG. 


(21) 


(See  the  discussion  in  Modest  16],  page  512.) 

For  the  boundary  condition  at  the  propellant  surface,  we  require  that  such  a  boundary  condition 
is  satisfied  in  an  integral  sense  (see  the  discussion  in  Modest  [16],  p.  470).  The  analysis  starts  by 
writing  for  the  intensity  along  the  propellant  surface: 


I(r  =  rWy  s)  =  Iw(rw,s),  (22) 

for  n  ■  s  >  0  (i.e.,  the  intensity  leaving  a  surface  must  be  prescribed  in  some  fashion  for  all  outgoing 
directions  n  s  >0).  This  relation  is,  of  course,  general  in  that  it  applies  in  all  cases.  Using 
assumption  4,  that  the  surface  emits  and  reflects  diffusely,  the  preceding  reduces  to 

I(r  =  rw,s)  =  Iw(s),  for  n  s  >  0.  (23) 

Evaluating  (14)  on  the  surface,  multiplying  by  s  •  n  and  integrating  over  the  hemisphere,  yields 


Let 


L 


n-s>0 


Iw{s)s  •  ndi 2 


1 

47 r 


/  (° 

Jji  a>0 


4-  3<7  •  s)s  •  ndi 2. 


q  s  =  qt\  sin  0  cos  0  4  qt2  sin  0  sin  0  4  qn  cos  6 , 


(24) 

(25) 
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where  (. )  are 
cos#.  Substitution 


the  tangent  components  and  qn  the  normal  component,  and  recall  that  s  n  — 
yields 


L 


T18>  0 


Iw(s)s  ■  ndtt  =  -( G  +  2 qn), 
4 


(26) 


or, 


G  4-  2 qn  —  4 


L 


Iw{s)s  •  ndQ. 


n  s>0 


(27) 


Substituting  assumption  4  for  the  intensity  in  terms  of  the  surface  radiosity  yields 


G  +  2  qn  —  4  Ju.  (28) 

In  addition  to  this  relation,  we  require  that  the  normal  radiative  heat  flux  is  the  sum  of  the  incident 
and  reflected  (radiosity)  contributions: 


Qn  —  J\v  H.  (29) 

The  incident  radiation  is  eliminated  in  favor  of  the  radiosity  using  (9)  and  the  hemispherical  reflec¬ 
tivity  is  related  to  the  absorptivity  a  by  the  Kirchoff  law  for  diffuse-gray  surfaces. 
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I 

II 

<Q. 

(30) 

to  yield 

enlb  —  aJw 

Qn  =  1 

1  —  a 

(31) 

We  use 

(28)  and  (31)  to  eliminate  Jw,  yielding 

2<ln  =  0  [467 r/h  oG] . 

(32) 

Substituting  (20)  to  eliminate  q  yields  the  proper  boundary  condition  along  the  propellant  surface 


-2n  •  VG 

(Ks  +  KaW-AxiV) 


1 

2  —  n 


[4e7r4  -  qG]. 


(33) 


This  is  an  implicit  boundary  condition  for  the  unknown  function  G  (i.e.,  a  Robin  boundary  condi¬ 
tion). 

The  other  boundary  conditions  assume  either  periodicity  for  randomly  packer  propellants  or 
zero  normal  fluxes  for  sandwich  propellants.  G  is  periodic  at 


x  =  ±LX ,  y  =  ±Ly, 


d  G  n  T  dG 

—  =  ()  at  x  =  ±LX,  — 
ox  ay 


0  at  y  =  ±Ly, 

and  far  awray  from  the  propellant  surface  the  radiation  field  is  in  equilibrium 


(34) 

(35) 


q  — ►  0,  i.e.,  VG  =  0  as  2  — >  oc.  (36) 

To  summarize,  we  solve  (19)  for  the  incident  radiation  G,  subject  to  the  boundary  conditions 
(33-36).  Once  G  has  been  determined,  the  radiative  heat  flux  q  can  be  found  from  (20).  The 
unknown  parameters  that  must  be  modeled  are  Ks ,  Ka,  Kem,  A\,  a,  and  e. 
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3.2  Governing  Equations 

Radiation  affects  the  combustion  field  through  the  radiative  heat  contribution,  V  q  =  —  V  qr,  and 
it  affects  the  connection  conditions  at  the  propellant  surface  because  of  the  term  qn  in  (31).  The 
propellant  surface  is  considered  to  be  an  opaque  surface,  in  the  sense  that  no  radiation  penetrates 
into  the  solid  phase  (i.e.,  all  is  absorbed  at  the  surface).  Thus,  no  changes  are  necessary  to  the  solid- 
phase  heat  conduction  equation.  We  note,  however,  that  the  presence  of  aluminum  does  affect  the 
heat  conduction  in  the  solid,  mainly  because  of  changes  in  thermal  conductivity  of  the  aluminum 
fuel  blend.  This  effect  is  accounted  for  by  using  homogenization  formulas  as  discussed  by  Jackson 
et  al.  [17].  These  formulas  will  be  described  in  more  detail  below. 


3.2.1  Gas  Phase 


In  this  section  we  briefly  outline  the  gas-phase  equations  to  highlight  the  changes  due  to  the  presence 
of  the  radiative  field. 

Consider  a  reacting  mixture  composed  of  Ns  species  that  can  react  according  to  Nr  reaction 
steps.  Let  T  denote  the  gas-phase  temperature  and  >*,  i  =  1  ,  ...,ATS,  the  gas-phase  species  mass 
fractions.  The  equations  expressing  conservation  of  energy  and  species  mass  in  a  non-conservative 
form  are 


DT  dP0 
pCp  Dt  ~  dt 


Nr 

~  V  q  4 


(37) 


p15t  =  v '  (^Wfc)  +  %  i3k'lRi’  k  =  1 . Ns'  (38) 

where  and  ai  axe  coefficients  based  upon  mass  and  total  enthalpy  conservation,  respectively, 
that  are  fully  determined  given  a  particular  reaction  scheme.  The  reaction  steps  are  expressed  in 
general  form  as 

Ri  =  DgyiP?Tm*  n  yf'fcexp(-^/T),  (39) 

k=  1 


where  Dq>t  is  the  gas-phase  Dainkohler  number  for  species  i;  Po  is  the  leading  order  pressure  in 
a  small  Mach  number  expansion  and  can,  at  most,  be  a  function  of  time;  and  Oi  the  gas-phase 
activation  temperature  for  species  i.  For  global  kinetics,  the  species  exponent,  in  the  reaction 
rates  need  not  be  related  to  the  stoichiometric  coefficients  of  the  corresponding  reaction;  their 
values  are  assigned  by  matching  reaction  rate  constants  with  experiments.  The  parameters  of  the 
combustion  model  {Dg^,rii,mi,vlk,  and  0{)  are  calibrated  using  a  global  optimization  strategy 
discussed  bv  Massa  et  al.  [18].  The  values  of  the  combustion  model  parameters  are  identical  to 
those  reported  in  Jackson  et  al.  [17]. 

In  addition,  p  is  the  gas-phase  density,  Cp  the  specific  heat  assumed  to  be  independent  of  the 
composition.  A  the  gas-phase  thermal  conductivity  taken  to  be  a  function  of  the  temperature  and 
the  volume  fraction  of  aluminum  and  aluminum  oxide  (AI2O3)  in  the  gas-phase,  and  q  the  heat-fiux 
vector.  In  writing  down  these  equations  the  Lewis  number  for  each  species  was  assumed  to  be  unity, 
and  second-order  or  higher  terms  in  the  Mach  number  were  neglected.  The  heat-fiux  vector  is  the 
sum  of  a  conductive  part  and  a  radiative  part. 


Q  =  Qc  +  Qr 


(40) 
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where  the  conductive  part  is  assumed  to  obey  Fourier’s  law 


qc  =  -A  VT, 


(41) 


and  the  radiative  part  is  given  in  the  previous  section 

-1 


Qr 


( Ks  +  Aa)(3  —  A\lo) 


VG, 


and 


v  qr  =  4nKemIb  -  KaG. 
Recall  that  If,  is  the  black  body  radiation 

h=-T\ 

7 r 


(42) 

(43) 

(44) 


and  the  incident  radiation  term  G  is  governed  by  the  equation  (19). 

The  pressure,  P,  is  written  as  the  sum  of  two  contributions,  a  leading-order  term  and  a  contri¬ 
bution  of  order  0(M2)  where  M  is  the  Mach  number, 


P  =  Po  +  P ,  tj-  =  0(M2).  (45) 

The  leading-order  term  Po  is  also  called  the  thermodynamic  pressure,  and  the  second  term  p  is 
called  the  hydrodynamic  pressure.  The  Navier-Stokes  momentum  equation  then  becomes 


VP0  =  0 


>  Po  =  Po(t)* 

Du  _ 
p— -  =  V  •  <7, 


(46) 


where  u  =  (w,  v,w)  are  the  velocity  components,  the  mean  pressure  Po  is  at  most  a  function  of 
time,  and  g_  is  the  total  stress  tensor,  the  sum  of  the  hydrodynamic  pressure  and  viscous  stress 
tensors, 

q_  =  -  ( p+  -/LtV  •  u  j  ]_  +  //  4-  (Vu)Tj  ,  (47) 


where  /  is  the  identity  tensor,  //  is  the  gas-phase  viscosity  and  is  related  to  the  gas-phase  thermal 
conductivity  via  //  =  Pr\/cv ,  and  the  superscript  T  denotes  the  matrix  transpose.  A  constant 
value  for  the  Prandtl  number,  Pr  =  0.72,  is  assumed. 

The  continuity  equation  is 

Dp 

'  •  +  pV  ■  u  =  0,  (48) 


and  the  equation  of  state  reads 


Po 


pRvT 
Mw  ’ 


(49) 


where  Pu  is  the  universal  gas  constant  and  Mw  is  the  averaged  molecular  weight  of  the  gas. 

Although  our  code  can  include  the  Navier-Stokes  equations,  in  this  work  we  use  the  Oseen 
approximation  which  significantly  reduces  the  computational  burden,  and  calculations  for  sandwich 
propellants  attest  to  its  accuracy  [19].  The  Oseen  approximation  sets  the  mass  flux  nominally 


11 


normal  to  the  burning  surface  to  a  constant  and  the  mass  flux  nominally  tangent  to  the  surface 
to  zero,  thus  satisfying  the  continuity  but  jettisoning  the  momentum  equation.  An  alternative 
and  essentially  equivalent  formulation  is  to  set  the  density  to  a  constant,  so  that  the  equation  of 
state,  Charles’s  law,  is  jettisoned;  and  a  uniform  velocity  field  u  =  0  and  v  =  constant  is  adopted, 
which  satisfies  both  the  continuity  equation  and  the  momentum  equation.  Then  the  mass  flux  pv 
is  constant,  and  the  model  is  merely  a  variation  on  that  proposed  by  Burke  and  Schumann  in  1928 
[20].  This  approximation  has  proven  to  be  adequate  and  suffices  for  our  purposes  here. 


3.2.2  Solid  Phase 


If  the  effects  of  mechanical  stress  are  neglected,  and  if  condensed  phase  reactions  are  neglected  in 
lieu  of  a  reaction  sheet  at  the  surface,  the  solid  phase  can  be  characterized  by  its  temperature  field. 
The  temperature  distribution,  Tc,  inside  the  solid  phase  is  governed  by 


PcCp- 


OTc 

dt 


=  V  •  (ACVTC), 


(50) 


where  pc  is  the  density  of  the  solid,  cv  the  specific  heat  assumed  to  be  equal  in  the  gas-phase,  and 
Ar  the  solid  phase  thermal  conductivity.  The  density  and  thermal  conductivity  are  functions  of  the 
spatial  coordinates  and  are  related  to  the  propellant  morphology  by  means  of  a  particle  marker 
method 

Pap  a  =i  ^AP  ^  =  1 

Pblend  \  A^lenci  ^  =  —  T 

where  ij)(x,y,z)  is  a  particle  marker  flag  such  that  a  point  (x,y,z)  lies  in  the  AP  if  ^  —  1,  and  in 
the  blend  if  ip  =  —1.  The  blend  can  be  a  homogenized  mix  of  aluminum,  AP,  and  binder.  The 
density  of  the  blend  is  then  given  by 


Pc 


(51) 


Pblend  =  ^APPAP  +  ^AlPAl  +  ^BPB,  (52) 

where  p\p  =  1.95  gm/cc  is  the  density  of  AP,  p\\  =  2.70  gm/cc  is  the  density  of  the  aluminum, 
pB  =  0-92  gm/cc  is  the  density  of  binder,  and  U  ( i  =  AP,  Al,  B)  is  the  volume  fraction  of  the  blend 
that  is  the  component  i.  Here,  t\p  +  £ai  +  ^B  =  1* 


3.2.3  Pyrolysis  Laws 

The  pyrolysis  laws  for  AP  and  binder  are  given  by  the  forms 

n,AP  =  -4ap  e~E^/RuTs ,  r(,.B  =  -4b  e~EB/R',Ts ,  (53) 

where  Ts  is  the  local  surface  temperature.  Since  aluminum  does  not  have  a  pyrolysis  law,  for  the 
blend  we  take 

rb, blend  =  -4 bien(]  e~Ebl^RuT’ ,  (54) 

where 

-4blend  —  (^Ap)^AP  (Afi)*8,  £blend  =  ^AP^AP  +  ^B^B^  (55) 

with  £ap  +  ^Al  +  =  1.  It  is  not  possible  at  this  time  to  verify  the  correctness  of  (54)  and  (55) 

since  a  fully  functional  three-dimensional  code  with  discrete  aluminum  is  not  complete,  and  so  we 
use  these  relations  with  some  reservations.  We  hope  to  return  to  this  point  at  a  later  date. 
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3.2.4  Propagation  of  the  Surface 

Suppose  the  surface  is  defined  by 

r}{x,y,z,t)  =  0.  (56) 

A  point  on  the  surface  moves  along  the  normal,  and  i]  does  not  change  for  such  a  point,  so  that 


Dri  ^  ( dx  dy  dz 

—  =0,  or  %  +  <-V,  =  0, 


The  normal,  pointing  into  the  gas,  is 


n  = 


Vr; 

Wv\' 


and  the  speed  of  the  surface  regression  rb  (>  0)  is  related  to  v  and  n  by 

v  =  - rbn ,  |v|  =  rb; 

rb  =  j'b(x,z,t)  is  the  local  burning  rate  of  the  solid.  Then  (57)  is 

fjt  -  rbn  •  Vr?  =  0, 
or 

Vt  ~  n\^v\  =  0. 

This  is  the  kinematic  equation  that  governs  the  propagation  of  the  surface. 

3.2.5  Connection  Conditions 

The  connection  conditions  at  the  propellant  surface  relate  the  solution  in  the  solid  phase  to  that 
in  the  gas  phase.  The  burning  interface  is  treated  as  a  reacting  regressing  sheet  separating  the  two 
phases.  The  temperature,  the  normal  mass  flux,  and  the  tangential  velocity  are  continuous  across 
the  sheet,  so  that  if  n  is  the  surface  normal  pointing  into  the  gas, 


(57) 

(58) 

(59) 

(60) 
(61) 


[T]  =  0, 

[M]  =  [p{u  ■  n  +  rfc)]  =  0. 

[«  x  n]  =  0. 

The  energy  balance  and  the  species  mass  balance  across  the  interface  become 

[-n  •  q\  -  —QSM, 


—n  ■  VYi 


=  M[Yi\. 


(62) 

(63) 

(64) 

(65) 

(66) 


The  jump  operator  in  the  above  equations  is  defined  as  [•]  =  (-)g  —  (-)c.  If  we  eliminate  the  heat  flux 
vector  in  the  energy  balance  in  terms  of  the  corresponding  conductive  and  radiative  parts,  then  we 
get 

[An  •  VT]  =  n-qr\g  -  QSM,  (67) 
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or,  upon  using  (32), 


[An  •  VT]  =  — -(4enlb  -  aG)  -  QSM . 


(68) 


We  note  that  the  phase-change  heat  release,  Qs,  is  a  function  of  space  along  the  surface  according 
to 


(  Q,s,AP  '0=1 

\  Qs, blend  tft  =  —1. 


(69) 


For  the  blend  we  have 

Qs,  blend  =  WAJpQs,AP  +  ^’AlQ,s,Al  +  UjbQs,B,  (70) 

where  Wi  ( i  —  AP,A1,B)  are  the  weight  fractions  of  the  blend  that  is  the  component  i .  Here, 
wap  4-  wai  +  w&  =  1.  Since  aluminum  does  not  pyrolyze,  we  set  Qs,a\  =  0. 


3.2.6  Coordinate  Mapping 

Roc  fire  uses  a  mapping  that  maps  the  non-fiat  moving  interface  onto  a  computational  stationary 
fiat  plane.  If  we  assume  the  surface  to  be  single- valued,  then  we  use 


v  =  y-  f{x,z,t). 

The  normal  and  tangent  vectors  at  the  surface  are  given  by 


n  = 


vi  +  /?  +  /!' 


t\  •  n  =  0,  t<2  =  n  x  t\, 


with  the  constraints 


|til  =  N  =  i. 

Derivatives  in  the  transformed  space  are  replaced  in  the  following  fashion 

o_  ^  o__  d_ 

dx  dx  x  Oij ' 

a_  d_ 

dy  dr] ' 

d  .  d 


0_ 

Oz 


dz  ^ 2  drj  ’ 


D_ 

Dt 


d_  ^  d_  _  d_ 

dt  dt  1  dy 

d  d  d  d 

m+Udi+Vfr1+WTz’ 


where 

and 


V-V-  ufx  -  wfz  -  ft, 

d2  /,  ,2  r2\  d2  &  „d2 

—  + (l  +  fx +  /,)-  +  —  - 2fx—~  2 fz 


d 


(71) 


(72) 

(73) 


(74) 


(75) 

(76) 
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Note  that  v  is  now  eliminated  as  a  dependent  variable  in  favor  of  v,  and  the  boundary  condition 
(64)  can  be  rewritten  as 


u  =  ~fx 


v  +  wfz  +  ft 


u>  =  ~fz 


1  +  /I  )'  ^  \  1+/? 


v  4-  ufx  -b  ft 


(77) 


In  the  Oseen  approximation  we  set  u  =  w  =  0  and  so  v  =  —  ft. 

3.2.7  Grid  Stretching 

It  is  convenient  to  first  reflect  the  solid  domain  about  the  77  —  0  plane  so  that  all  equations  are 
defined  on  the  domain  r/  >  0.  Once  this  has  been  done,  an  additional  grid  stretching  mapping  can 
be  used  to  cluster  the  computational  points  in  the  //-direction  near  the  surface.  The  grid  clustering 
mapping  is  given  by  the  functional  relationship  fj  —  77(77).  Then 


and 


dr]2 

Roc  fire  currently  uses  the  mapping 


dr]  0 
dr]  dfj  ’ 

(78) 

2  d2  d2fj  d 
dfj 2  +  dr)2  df) 

(79) 

V 

WriPf’ 

(80) 

7}  = 


where  Lv  is  the  end  of  the  computational  domain  in  the  77-direction,  and  C  is  a  clustering  constant  . 
Typically  we  choose  C  —  6. 

3.2.8  Transformed  Equations 

In  this  section  we  summarize  the  equations  and  boundary  conditions  in  the  transformed  (x,7],z) 
space,  with  0  <  fj  <  +00. 

77  >  0 


For  the  temperature,  species,  and  radiation  equation  we  have 
dfj  \  dP{ 

+ wTf  ~  Tt 

( 

di] 


pcv  (  Tt  +  uTr  +  v  ^Tn  +  wTz)  -  ^  =  V  •  (A  VT)  -  (4nKemIb  -  KaG )  +  jP  a{Ru 

1=1 


P  (  YkJl  +  uYk,x  +  f,  f  yM  +  wYk,z)  =  V  •  ( - vrfc )  +  fajRi ,  k  =  1 . Ns, 


Nr 

Z 

1  =  1 


V  di]z  \  dr) 


GXx  +  (1  +  fx  +  fz)  (  +  I  ~r  )  Gfjfl  I  +  Gzz 


dr] 


(81) 

(82) 

(83) 


^  i^fxGxt]  +  2 fzGZjj  +  (fx x  +  fzz)Gff)  —  (As  +  Aa)(3  —  A\oj)  RqG  —  . 
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where 


v  •  (wt) = (i + a + ih  (0(Ari) + (^)5  mu) 

+t\T,U  +  (AT;)*  -  ^  </,!AT,),  +  /,(AT,),  +  /,(AT„),  +  MAT,),) , 

and  similarly  for  the  species  Y*. 

For  the  fluid  equations,  they  remain  unchanged  due  to  radiation,  so  they  are  not  repeated  here. 


For  the  condensed  phase  temperature  we  have 

PcCp  (tc4  +  ft%±  Tc^  =  (1  +  fl  +  ft)  (^(X  cTc,n)  +  (^)2  (A  cTc,n)^ 

“I"(^c7c,x)x  4-  (A rTc,z)z  +  (/x(AfTc  ^)x  4-  /x(A cTc,x)r)  4-  /z(AcTc^)2  4-  /2(ArTr  2)^) . 
Note  that  (85)  is  the  equation  used  in  our  Rocburn  code. 


(85) 


T)  — >  OO 

4(.)  =  0,  Tc  =  To,  (86) 

where  To  is  the  cold  supply  temperature  deep  in  the  solid  phase,  and  the  notation  (•)  means  the 
quantities  (Tg,Tc,YuiLV,w,p). 


Tj  —  0 


u  =  -/: 


/t+Wl  +  /|  +  /?, 
r(x,(),2,0lg  =  T(x,0,2,f)lc. 
pgv  =  pcrb  V  1  +  Si  +  fl  =  M  V7!  +  /x  +  /f, 

fl  +  wfz  +  ft\  _  t  A  +  nfx  +  ft  \ 


1  +  /J  )' 


tc 


Aff  (  -/,r*  +  (1  +  /|  +  fl)~^Tn  -  fzTz 


1+/?  ; 

+  Ac  ( /,r,  +  (1  +  fl  +  /22)$ r,  +  /2r2 


dr? 


1 


2(2  -  a) 
dpy 
*'dn  *'n 


-fzYi*  +  (1  +  fl  +  fl)T-YiA  -  fzYi 


(4«r4  -  ftG)  -  Q,Mj  ^1  +  ft  +  ft, 

=  m  (r,|„  -  v. 


-2(2 -ft) 


(As  +  Aa)(3  —  .Aia;) 


-/xGx  +  (1  +  fl  +  ~  fzG 

=  (467 Th  -  oG) \g  v/1  +  Tl  +  ft. 


where  K  ,  denote  fix  surface  values. 


(87) 

(88) 

(89) 

(90) 

(91) 

(92) 

(93) 
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•E  —  Dxt  Z  —  T  Lz 


Periodicity. 


3.3  Radiation  Parameters 


The  radiation  model  is  an  adaption  of  the  model  used  by  Brewster  and  Parry  [15]  for  conventional 
size  aluminum,  and  used  in  our  previous  study  of  a  blend  of  fine  AP  and  aluminum  [17].  Thus  the 
contribution  of  the  particles  is  accounted  for  by  setting, 


oxide  &  oxide  >  ^Ai^Al 

I<a  =  — f; - +  *  1 


K»  = 


KP 


D oxide 
♦^oxidePoxide 

•Doxide 
dioxide  f  oxide 
^oxide 


k\ 


+  k\ 


Dm 

df  A1PA1 

Dm  ' 

df  AKAl 

Dm  ' 


(94) 


where  tM  and  ^oxide  are  the  volume  fractions  of  aluminum  and  aluminum  oxide,  respectively,  in 
the  gas-phase;  a,  p,  and  e  are  the  absorptivity,  reflectivity,  and  emissivity,  respectively;  and  k\  is 
an  empirical  constant  that  reflects  uncertainties  in  the  evaluation  of  the  diameter  of  the  aluminum 
particles.  In  the  simulations  discussed  in  this  paper,  this  parameter  has  been  set  to  unity. 

The  average  diameter  of  the  injected  particles  can  be  much  larger  than  the  average  diameter 
of  the  initial  distribution  of  aluminum  due  to  the  possibility  of  agglomeration.  Agglomeration  is  a 
phenomenon  that  depends  on  the  aluminum  particles  surface  residence  time,  which  in  turn  depends 
on  the  propellant  burn  rate  [21].  For  burn  rates  typical  of  composite  propellant  with  ultrafine 
aluminum,  the  residence  time  is  short  and  agglomeration  effects  must  be  considered  negligible. 
The  experimental  investigation  of  Stephens  et  ai  [5]  has  shown  that  changing  the  size  of  aluminum 
from  3  microns  to  nano-size  aluminum  does  not  alter  the  burn  rate  for  monomodal  200  /mi  diameter 
AP.  However,  the  addition  of  36-//m  aluminum  sharply  reduces  the  propellant  burn  rate  when  the 
overall  volume  fraction  of  metal  is  maintained  constant.  Consequently,  all  simulations  discussed 
in  this  work  assume  that  the  initial  diameter  of  the  particle  ejected  in  the  gas  phase  is  Do  = 
3-/mi  corresponding  to  the  upper  boundary  of  the  fine/ ultrafine  range  as  it  was  defined  in  the 
introduction. 

Other  model  uncertainties  include  the  values  of  oai,  Pa\,  and  6ao  as  well  as  the  diameter 
D0x ide-  The  values  of  the  constants  used  in  this  study  are  given  in  table  (3),  and  are  taken  from 
Brewster  and  Parry  [15].  The  diameter  of  aluminum  and  oxide  particles  are  determined  by  solving 
for  Eulerian  fields  representative  of  the  dispersed  gas-particulate  mixture;  this  topic  is  discussed  in 
detail  in  the  next  section. 

The  parameters  listed  in  equation  (94)  are  identical  to  those  used  by  Jackson  et  ai  17  and 
consider  the  contribution  of  the  particles  only.  The  contribution  of  the  gas  mixture  is  accounted 
for  by  evaluating  the  effective  absorption  coefficient  for  the  gray  gas  as  the  Planck-mean  absorption 
coefficient,  Kv,qaS)  Modest  [16],  page  344.  Kpms  should  be  in  principle  determined  by  integrating 
the  spectral  absorption  coefficient  of  the  gas  over  the  relevant  bands.  The  result  of  this  calculation 
would  be  temperature,  pressure  and  composition  dependent.  Due  to  the  absence  of  a  detailed 
kinetic  mechanism  in  our  model,  only  the  pressure  dependence  is  retained, 


Dp, gas  —  ftp, gas  D. 


(95) 
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Assuming  the  combustion  products  to  be  composed  in  large  part  of  CO 2,  we  select  the  value  Kvms  — 
0.3  cm  atm-1,  based  upon  the  results  of  Zhang  and  Modest  22] .  For  the  gas-particulate  mixture, 
the  coefficients  are  determined  by  augmenting  Ka  and  Kem  in  equation  (94)  by  the  contribution 
hp,gas- 


a 

P 

e 

D  (/zm) 

Al 

0.10 

0.90 

1.00 

1.0 

AI2O3 

0.45 

0.55 

0.45 

0.3 

Table  3:  Radiation  parameter  values;  from  Brewster  and  Parry  [15]. 


3.4  Aluminum  Modeling 

The  presence  of  aluminum  in  the  solid  phase  is  accounted  for  by  evaluating  the  density,  ther¬ 
mal  conductivity,  and  pyrolysis  rate  according  to  homogenization  formula  derived  for  line  binder- 
aluminum  blend.  Most  of  the  homogenization  formulae  are  dictated  by  conservation  principles, 
and  their  derivation  is  straightforward.  The  expression  for  the  thermal  conductivity  does  not  stem 
form  a  conservation  argument.  The  formulae  used  in  this  research  are  identical  to  those  described 
by  Jackson  et  al.  [17], 


n3  ( _ 1  (Aai/AhTPB) 

Al  V  (-^blend/'^HTI’B)  ~  (Aai/AhtPb) 


3 

(Ablend/AHTPB) 


(96) 


where  blend  identities  the  fuel-aluminum  mix  in  the  solid.  Likewise,  in  the  gas  phase,  the  thermal 
conductivity  is  evaluated  from  homogenization  formulae,  which  are  formally  identical  to  those  in  the 
solid  phase,  equat  ion  (96),  and  require  knowledge  of  the  gas  and  particulate  thermal  conductivities 
and  gas  phase  volume  fractions. 

The  presence  of  aluminum  in  the  gas  phase  is  modeled  by  defining  a  set  of  continuous  Eulerian 
fields  that  describe  the  particle  transport  as  a  sub-scale  phenomenon  of  the  gas  phase  micron-scale 
processes:  those  based  upon  morphological  length  scales. 

The  aluminum  in  the  gas  phase  is  present  either  as  liquid  aluminum  or  solid/liquid  aluminum 
oxide.  These  two  states  are  connected  along  the  global  chemical  path, 

Al(s)  +  mOX  — *  ^Al203(s//)  +  products,  (97) 


where  OX  typifies  oxidizing  species  in  the  product  stream  and  m  is  equal  to  3/2  divided  by  the 
number  of  oxygen  atoms  in  OX.  The  aluminum  reacts  with  the  hot  oxidizing  species  away  from  the 
surface  because  of  the  inhibiting  effect  of  the  oxide  shell.  Therefore,  OX  includes  the  oxygen  rich 
species  in  the  equilibrium  gas  products  of  the  AP-HTPB  combustion.  If  we  assume  the  specific  heat 
of  reactants  and  products  equal,  the  heat  release  can  be  found  based  upon  the  heat  of  formation 
of  the  species.  Hence,  if  OX  =  O2  in  equation  (97),  the  associated  heat  release  is  Q^\  =  7165 cal/g, 
while  if  OX  =  H2O  the  heat  release  is  Qai  —  3957 cal/g\  note  that  this  value  is  per  gram  of 
aluminum  and  that  the  aluminum  oxide  is  assumed  in  the  liquid  state.  Based  upon  equilibrium 
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chemistry  calculation,  by  matching  flame  temperatures  for  a  AP-HTPB-A1  homogenized  blend,  we 
select  the  value  to  be  closer  to  the  latter,  and  we  set  QA\  =  4500ca//g. 

Given  that  the  fields  are  modeled  as  Eulerian  fields,  we  introduce  at  this  point  the  concept 
of  Eulerian  density  of  the  nano-scale  condensed  phases,  and  distinguish  it  from  the  Lagrangian 
density,  which  is  constant  and  equal  to  the  material  value.  The  Eulerian  density  is  defined  as  the 
ratio  of  the  mass  of  condensed  matter  to  the  extension  of  the  volume  that  contains  it,  when  the 
volume  is  taken  smaller  than  all  the  relevant  physical  scales.  The  Eulerian  densities,  pA\,  Poxide  are 
related  to  the  Lagrangian  counterparts,  PAl.oxide  by, 


PAI, oxide  —  P  Al,  oxide  ^Al,  oxide? 


(98) 


where  £ai, oxide  are  the  volume  fractions  of  aluminum  and  aluminum  oxide,  respectively,  in  the  gas 
phase.  Also  we  note  that  the  overall  density  is, 


P  —  Poxide  “h  PAI  Pgas? 

where  the  density  of  the  gas  is  found  from  the  Dalton’s  Law  for  a  mixture  of  perfect  gases, 

P 

Pgas  (1  —  ^A1  —  f  oxide)  • 

Finally,  the  mass  fractions  of  the  dispersed  particle  phase  are, 


Yi 


A 1, oxide 


PA1, oxide 


(99) 


(100) 


(101) 


The  preceding  formulas  yield  a  modified  thermal  equation  of  state  in  which  the  mass  fractions 
of  the  particulate  field  modulate  the  pure  gas  equation.  Specifically,  in  terms  of  the  Lagrangian 
densities  of  the  particulate,  its  mass  fractions,  the  chamber  pressure  and  the  temperature. 


p  Rt{1  +  Ym  (  RTpM  0  +  r°xide  V  RTpox ide 
The  transport  equations  for  the  Eulerian  fields  are 


-l 


DYa  i 


oxide 


Dt 


^  *  (^Al,oxide^^  Al, oxide)  +  ^Al, oxide* 


(102) 


(103) 


where  V  is  a  diffusion  coefficient  which  models  the  motion  of  the  particles  with  respect  to  the  gas 
flow  due  to  concentration  gradients,  a  phenomenon  necessary  to  maintain  spatial  continuity  of  the 
concentration  fields.  Its  value  is  assumed  to  be  very  small,  and  taken  to  be  the  fifteenth  part  of  the 
diffusion  coefficient  of  the  gaseous  species.  This  is  a  convenient  value  because  it  corresponds  to  the 
minimum  value  necessary  to  obtain  convergence  of  the  Mult  igrid  Poisson  solver  for  the  radiative  heat 
equation.  Simulations  have  demonstrated,  however,  insensitivity  to  this  parameter  and  therefore 
it  has  not  been  selected  for  a  parametric  study.  The  term  uj  in  equation  (103)  represents  the  mass 
conversion  of  aluminum  to  alumina  and  is  determined  from  the  aluminum  burning  model  discussed 
in  detail  in  the  following  section.  Note  that,  once  the  mass  fractions  are  calculated,  the  diameter 
of  the  particles  is  evaluated  as, 


DX\  =  D0 


f  YAlflM  \ 

~h  4  oxidp  / 


1/3 


(104a) 
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1/3 


^oxide  —  Do 


_ ^  oxide _ 

^  A\  ft M  “I"  ^oxide  ftp 


(104b) 


where  ftp  =  p0xide/PAb  while  ftw  =  Mox\c\e/(2MA\ ),  with  M  the  molecular  weight.  D0  is  the  initial 
diameter  of  the  particles  as  they  are  ejected  from  the  surface. 


3.5  Aluminum  Burning 
3.5.1  Model 


The  burning  of  aluminum  spheres  is  modeled  using  a  power  law  representation  of  the  volume 
consumption  of  the  particles  in  a  constant  property  environment, 


Ov 


l_dV 
V  dt 


-k  D~n 


(105) 


where  the  exponent  n  is  commonly  assumed  to  be  close  to  2  and  the  cofactor,  kv,  is  dependent 
on  the  local  conditions,  such  as  temperature,  pressure  and  oxidizing  species  concentration.  The 
exponent  and  cofactor  can  be  obtained  from  experimentally  measured  particle  burn  out  times  if 
the  surrounding  conditions  are  assumed  invariant  of  time.  Experimental  correlations  are  usually 
presented  as 

turnout  =  kbDgdT~n‘ P~n” ,  (106) 

where  xox  is  the  concentration  of  oxidizing  species  active  along  the  path  equation  (97)  and  VTe\  is 
a  diffusion  parameter,  see  the  discussion  in  Najjar  et  al  [23]  for  more  details.  Time  integrating 
between  the  initial  time  and  the  burn-out  time  yields  the  relations. 


n^nd  (107) 

K  =  ±-r*«pn  p^Drel.  (108) 

The  value  of  obtained  by  different  experimental  procedures  can  differ  by  several  orders  of  mag¬ 
nitude.  No  power  law  exists  that  is  valid  in  all  the  range  of  aluminum  diameters  of  interest  to  solid 
propellant  analysis. 

Another  modeling  uncertainty  is  rooted  in  the  fact  that  aluminum  near  the  surface  burns  dif¬ 
ferently  before  and  after  the  oxide  shell  break-up.  For  micron-scale  particles,  the  shell  break-up 
is  commonly  identified  with  the  melt  temperature  of  the  aluminum  oxide,  Tmeit tOXide  —  2350 A";  for 
nano-scale  particles  the  formation  and  re-sealing  of  cracks  makes  this  definition  more  ambiguous. 
Below  such  transitional  temperature  the  micron-scale  aluminum  is  assumed  to  undergo  negligible 
oxidation.  On  the  contrary,  nano-scale  particles  are  considered  to  undergo  significant  oxidation  at 
temperature  below  that  of  the  shell  break-up,  at  approximately  the  aluminum  melting  temperature 
24,25].  Experiments  by  Bazyn  et  al  [26],  have  shown  that  the  burn-out  time  of  nano-scale  particles 
is  essentially  temperature  insensitive  for  ambient  value  lower  than  approximately  2000  K,  while  it 
sharply  decreases  of  an  order  in  magnitude  above  such  a  temperature,  demonstrating  the  presence 
of  two  stages  of  the  burning  in  nano-scale  aluminum.  The  first  stage,  for  T  <  Xineitioxide^  is  char¬ 
acterized  as  the  core  shrinking  phase  of  the  burning,  a  strictly  diffusion  limited  and  temperature 
independent  process,  Aita  et  al  [27];  the  reason  why  the  core  shrinking  phase  is  not  important  at 
the  micron-scale  level  is  the  larger  volume/surface  ratio  of  the  particles.  The  second  stage,  valid 
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after  the  oxide  shell  brakes,  is  more  complex  and  its  model  should  take  into  account  aluminum 
kinetics  paths. 

Widener  and  Beckstead  [28]  model  offers  the  best  available  correlation  for  the  full  combustion 
stage  of  aluminum  burning.  This  correlation  was  obtained  for  micron-size  aluminum  so  it  does  not 
necessarily  extend  to  the  nano-scale  range.  By  using  Widener  and  Beckstead  correlation  we  set, 
kb  =  1138,  nt  =  1.57,  np  =  0.2,  nQ  =  0.39,  when  the  dimensions  are  taken  to  be  degree  Kelvin, 
atmosphere  and  microns.  Due  to  the  significant  disparity  in  the  burn-out  times  between  the  two 
stages  of  burning,  the  core  shrinking  stage  is  disregarded,  by  effectively  setting  kv  =  0  when  the 
temperature  is  lower  than  the  melting  temperature  of  the  oxide. 

Experiments  by  Bazyn  et  al.  [26]  have  shown  that  the  temperature  overshoot  above  ambient  of 
aluminum  particle  in  constant  ambient  experiments  is  limited  by  the  dissociation  of  alumina.  To 
reproduce  this  phenomenon,  aluminum  burning  is  turned  off  when  the  gas  temperature  reaches  the 
boiling  point  of  the  AI2O3.  The  average  value  for  the  aluminum  oxide  boiling  point  in  the  pressure 
range  P  G  20-70  atm  is  taken  to  be  Tboii, oxide  —  4100 K.  These  considerations  yield  a  transformed 
particle  burning  law, 


(Tv 


—  k  D~n 


Tmelt, oxide  <  T  <  ^boil, oxide 
P  ^  ^rnelt, oxide  1  P  ^  Tboiroxide- 


(109) 


Equation  (109)  is  then  used  to  evaluate  the  source  terms  in  (103),  resulting  in  the  formalisms, 


^ai  =  pYk\Vv 


(110a) 


^oxide  —  ^  h\$\i  • 


(110b) 


3.5.2  Uncertainty  and  parametric  study 

Due  to  the  high  degree  of  uncertainty  in  the  burning  formulas,  a  parametric  study  will  be  conducted 
by  pre- multiplying  kv  by  a  factor,  denoted  in  the  following  as 

K  =  ^rn‘Pn'x£Prel-  (HI) 

The  factor  will  be  varied  by  effectively  calibrating  the  model  versus  experimental  data.  The  cali¬ 
bration  shows  <pv  to  be  always  less  than  one  in  order  to  obtain  sensible  results,  demonstrating  that 
micro-scale  burn  rate  correlations  over-predicts  the  burn  rate  of  nano-scale  particles. 


4  Results 

4.1  Validation  Study 

For  validation  we  consider  propellant  POla  of  [2].  The  formulation  for  this  propellant  is  15%  by 
weight  nano-aluminum,  17%  HTFB  binder,  and  68%  AP.  The  AP  is  biinodal  with  80%  by  weight 
150  micron  AP  and  20%  75  micron  AP.  The  morphology  is  modeled  using  Rocpack ,  our  proprietary 
packing  code.  A  3,124  particle  pack  that  models  POla  is  shown  in  Figure  (1).  Figure  (2)  shows  a 
grid  resolution  study  for  three  different  grids.  What  is  plotted  is  the  surface- aver  aged  burn  rate  as  a 
function  of  distance  burned,  and  for  three  different  grids.  Note  that  the  surface-averaged  burn  rates 
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are  not  constant  due  to  morphology  burn-through.  We  comment  here  that  the  grid  is  stretched  in 
the  nominal  ^/-direction,  so  only  60  grid  points  are  needed  to  adequately  resolve  the  combustion 
dynamics  as  well  as  the  far  field  conditions.  The  grids  in  the  (a\  ^-directions  are  uniform.  The  burn 
rates  for  the  three  grids  are:  0.9097  (40  x  40  x  60);  0.9649  (80  x  80  x  60);  1.0152  (160  x  160  x  60). 
It  is  seen  that  the  middle  grid  gives  adequate  resolution  to  within  5%,  and  so  we  use  this  grid  in  all 
subsequent  simulations.  Figure  (3)  shows  the  surface-averaged  burn  rate  as  a  function  of  distance 
burnt  for  various  values  of  the  uncertainty  parameter  (fiv  and  for  a  pressure  of  20  atm.  The  burn 
rates  are:  0.8568  cm/s  (0r  =  0.15),  0.8305  cm/s  (<; pv  —  0.08),  0.8137  cm/s  (cf>v  =  0.06),  0.7796  cm/s 
(<j>v  =  0.04).  The  value  for  (ftv  —  0.04  is  reasonable  close  to  the  value  reported  in  [2],  and  so  we  fix 
this  parameter  and  use  it  in  all  subsequent  simulations. 

Figure  (4)  plots  the  surface-averaged  burn  rate  as  a  function  of  distance  burnt  and  for  different 
pressures.  Note  that  the  surface- averaged  burn  rates  are  not  constant  due  to  morphology  burn- 
through.  Figure  (5)  compares  the  burn  rate  from  the  simulations  to  those  of  the  experimental  data 
of  [2].  For  both  the  simulations  and  the  experiments  we  also  show  the  error  bars.  Note  the  overall 
good  agreement  between  simulations  and  experiments.  Figures  (6-8)  plots  the  surface  temperature 
and  the  radiative  transport,  —V  •  qr,  for  three  different  times,  respectively. 


Figure  1:  A  3,124  particle  pack  that  models  POla  of  [2].  The  coarse  AP  particle  distribution  is 
shown  in  red,  and  the  fine  AP  particle  distribution  is  shown  in  gray.  The  void  space  is  the  blend, 
a  homogeneous  mixture  of  very  fine  AP,  binder,  and  nano-aluminum  particles.  The  domain  is 
periodic,  and  only  those  spheres  inside  the  cube  are  plotted.  The  cube  dimension  is  1079.79  x 
1079.79  x  1079.79  /zm3. 
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Figure  2:  Grid  resolution  study  showing  surface  averaged  burn  rates  as  a  function  of  distance  burnt 
at  a  pressure  of  40  atm  and  with  <j)v  =  0.04.  The  average  burn  rates  for  the  three  grids  are:  0.9097 
(40  x  40  x  60);  0.9649  (80  x  80  x  60);  1.0152  (160  x  160  x  60).  In  each  case  the  simulations  were 
run  on  40  processors.  The  dash  lines  correspond  to  the  length  of  the  periodic  pack.  Pack  shown  in 
Figure  (1). 
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Figure  3:  Surface-averaged  burn  rate  as  a  function  of  distance  burnt  for  various  values  of  <t>v  and  at  a 
pressure  of  20  atm.  The  average  burn  rates  are:  0.8568  cm/s  (<j>v  —  0.15),  0.8305  cm/s  (< pv  =  0.08), 
0.8137  cm/s  (4>v  =  0.06),  0.7796  cm/s  (<j)v  =  0.04).  The  dash  lines  correspond  to  the  length  of  the 
periodic  pack. 


Figure  4:  Surface-averaged  burn  rate  as  a  function  of  distance  burnt  for  various  pressures  (20,  40 
and  80  atm)  and  for  <j>v  =  0.04.  The  average  burn  rates  are  0.7796  cm/s  (20  atm),  0.9649  cm/s  (40 
atm),  and  1.2256  cm/s  (80  atm).  The  dash  lines  correspond  to  the  length  of  the  periodic  pack. 
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Figure  5:  Burn  rate  as  a  function  of  pressure  for  (pv  =  0.04.  Numerical  values  are  shown  as  red 
boxes  with  5%  error  bars.  Also  shown  are  the  experimental  values  from  [2],  shown  as  black  circles 
with  experimental  error  bars. 
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Figure  6:  Plot  of  surface  temperature  and  radiative  transport  at  time  0.090964  s. 
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(a)  Surface  Temperature 


(b)  Radiative  Transport 


Figure  7:  Plot  of  surface  temperature  and  radiative  transport  at  time  0.15296  s. 


(a)  Surface  Temperature  (b)  Radiative  Transport 

Figure  8:  Plot  of  surface  temperature  and  radiative  transport  at  time  0.18596  s. 
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4.2  Purdue  Propellants 

In  a  preliminary  study,  propellants  containing  AP  and  nanoaluminum  were  mixed  at  the  mixing 
facility  at  Purdue  University.  The  propellant  mixer  is  shown  in  Figure  (9).  The  baseline  batch, 
which  we  will  call  series  A,  is  based  on  historic  Purdue  formulations  of  micron-aluminum  loaded 
AP/HTPB  composites  13  with  a  direct  substitution  of  nano  aluminum  for  the  micrometric  alu¬ 
minum.  The  nanoaluminum  used  was  Novacentrix  80nm  aluminum  powder  (it  should  be  noted 
that  the  size  distribution  for  this  specific  batch  is  outside  of  manufacturer  specified  tolerances), 
and  all  other  materials  were  from  Firefox  Enterprises,  a  pyrotechnic  supply  for  amateur  hobbyists. 
Table  4  lists  the  mix  formulation. 


Table  4:  Series  A  Propellant  Formulation 


Ingredient 

Function 

Percent  Weight 

Hydroxyl-terminated  Polybutadiene  (R-45M) 
Ammonium  Perchlorate  (60-130rmnn  Diameter) 
Ammonium  Perchlorate  (200rmmi  Diameter) 

2  Ethyl  Hexyl  Acetate 

Isonate  143-L  -  Modified  MDI 

HX-878  -  Tepanol 

Lecithin 

Polydimethylsiloxane 

Binder  and  Fuel 

Oxidizer 

Oxidizer 

Plasticizer 

Curing  Agent 
Bonding  Agent /Curative 
Processing  Aid 
Anti-Foam  Agent 

14.84 

38.64 

25.00 

4.58 

2.95 

0.60 

0.50 

2  drops 

To  determine  the  effect  of  nano  aluminum  loading,  we  construct  three  propellants  with  the 
same  coarse-tofine  (C/F)  ratio  of  AP  as  the  baseline  batch,  but  vary  the  weight  of  nano-aluminum 
from  6%,  12%,  and  18%.  Table  5  shows  the  percent  weight  of  the  various  ingredients  used  in  the 
simulation  matrix.  Representative  packs  of  2,000  particles  are  shown  in  Figure  (10).  The  length 
of  each  side  of  the  cubes  are  given  in  Table  5.  The  surface-averaged  burn  rate  as  a  function  of 
distance  burnt  is  shown  in  Figure  (11).  Note  that  the  surface-averaged  burn  rates  are  not  constant 
due  to  morphology  burn-through.  The  burn  rate  for  the  lowest  aluminum  content,  Figure  (11a),  has 
larger  variations  about  the  mean  than  the  other  two  formulations,  which  indicates  a  more  erratic 
burning  behavior.  The  corresponding  burn  rates  are  shown  in  Figure  (12),  at  a  pressure  of  68  atm. 
The  black  squares  correspond  to  the  simulations  for  the  three  propellant  formulations  Al,  A2  and 
A3,  and  the  blue  circle  corresponds  to  the  experiments  carried  out  at  Purdue  for  propellant  mix 
A2.  The  95%  confidence  level  is  also  shown  in  blue.  The  burn  rate  values  from  the  simulations 
are  1.0725  cm/s  (Al;  6%),  1.1439  cm/s  (A2;  12%),  and  1.0902  cm/s  (A3;  18%).  Note  that  the 
largest  value  of  the  burn  rate  corresponds  to  12%  nano-aluminum.  The  experimental  burn  rate  for 
propellant  A2  (shown  as  the  blue  circle)  is  1.26  cm/s  (see  discussion  below).  This  corresponds  to  a 
10%  difference  between  the  simulation  and  the  experimental  burn  rate;  considering  the  error  bars 
are  probably  about  5%  for  both  the  experiments  and  the  simulations,  this  represents  an  excellent 
result. 

Due  to  the  small  amount  of  propellant  needed,  samples  were  prepared  by  hand-mixing  the 
chemicals  in  a  plastic  bag  as  opposed  to  the  Ross  mixer.  Despite  concerns  about  adequate  rheology 
and  particle  dispersion,  this  procedure  has  been  shown  to  provide  inconsequential  differences  6. 
Propellants  were  also  cast  in  an  oven  at  140°F  for  a  minimum  of  two  days,  and  cast  in  either 
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Mix 

nano-Al 

HTPB 

Coarse  AP 

Fine  AP 

L  (/an) 

A1 

6 

23.47 

27.734 

42.796 

1204.5 

A2 

12 

23.47 

25.374 

39.156 

1235.2 

A3 

18 

23.47 

23.015 

35.515 

1270.2 

Table  5:  Simulation  matrix  varying  percent  weight  of  nano-aluminum  with  fixed  C/F  AP  ratio  of 
0.6480331.  The  last  column  corresponds  to  the  length  of  one  side  of  a  periodic  cube  using  2,000 
total  particles. 


Figure  9:  Propellant  mixer  at  Purdue  University. 


(a)  6%  nano-al  (b)  12%  nano-al  (c)  18%  nanoal 

Figure  10:  Particle  packs  that  model  Series  A  from  Purdue.  The  coarse  AP  particle  distribution 
is  shown  in  red,  and  the  fine  AP  particle  distribution  is  shown  in  gray.  The  void  space  is  the 
blend,  a  homogeneous  mixture  of  very  fine  AP,  binder,  and  nano- aluminum  particles.  The  domain 
is  periodic,  and  only  those  spheres  inside  the  cube  are  plotted. 
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Figure  11:  Surface- aver  aged  burn  rates  as  a  function  of  distance  burnt  at  a  pressure  of  68  atm. 
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Figure  12:  Burn  rate  as  a  function  of  percent  nano-aluminum.  Pressure  is  68  atm.  The  black 
squares  correspond  to  the  simulations,  and  the  blue  circle  corresponds  to  the  experiments  carried 
out  at  Purdue.  The  95%  confidence  level  is  also  shown  in  blue. 
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Figure  13:  Propellant  cast  in  quartz  tubes. 


Figure  14:  Propellant  cut  from  cast  block. 


quartz  tubes  or  as  a  solid  block  and  cut  into  individual  strands,  as  shown  in  Figures  (13)  and  (14), 
respectively. 

Furthermore  the  cut  strands  are  inhibited  by  either  using  a  polyurethane  or  epoxy  coating  so 
that  the  propellant  can  only  burn  on  one  surface.  Figure  (15)  shows  a  fully  prepared  propellant 
strand  ready  to  go  into  the  6000  psi  Crawford  bomb. 

Initial  burns  with  the  Al,  A2,  and  A3  mixes  cast  in  quartz  tubes  yielded  uneven,  noil-planar 
burns  due  to  poor  packing  and  air  pockets  within  the  sample,  illustrated  in  Figure  (16).  Conse¬ 
quently,  burn  rate  analysis  is  erratic  and  inconsistent. 

A  redo  of  mix  A2,  cast  as  a  block  and  inhibited  with  epoxy,  resulted  in  much  more  planar 
burns,  shown  in  Figure  (17),  and  allows  for  consistent  surface  tracking  using  the  open  source 
program  Tracker  [33].  Analysis  of  four  strands  burned  at  an  initial  temperature  of  60.7°F  and 
initial  pressure  of  68  atm  resulted  in  a  burn  rate  of  1.26  cm/s. 
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Figure  15:  Inhibited  propellant  prepped  to  go  into  Crawford  bomb.  Shown  also  is  NiChroine 
ignition  wire. 


Figure  16:  Video  images  from  an  uneven,  nonplanar  burn. 


Figure  17:  Video  images  from  an  even,  planar  burn. 


5  Purdue  Experimental  Results 

In  efforts  to  firmly  validate  the  simulation  codes  developed  by  IllinoisRocstar  LLC,  a  test  matrix  of 
propellants  has  been  developed  at  Purdue  facilities  for  computer  model  verification.  The  following 
sections  detail  the  characterization  of  these  propellants,  the  experimental  approach  to  meet  the 
technical  objectives  for  this  project,  and  the  current  results  of  work  done. 

5.1  Propellant  Characterization 

For  the  simulation  codes  to  accurately  model  a  given  propellant,  information  must  be  known  about 
its  composition  and  particle  size  distributions  of  relevant  constituents.  Since  propellant  mass  com¬ 
position  is  easily  controlled,  a  simple  test  matrix  of  propellants  was  established  to  provide  validation 
and  calibration  for  the  IllinoisRocstar  code.  The  first  set  of  propellants,  designated  Series  A,  con¬ 
sists  of  nano-aluminum  (nAl)  fuel,  fine  ammonium  perchlorate  (fAP)  oxidizer,  coarse  ammonium 
perchlorate  (cAP)  oxidizer,  hydroxyl-terminated  polybutadiene  (HTPB)  binder,  and  various  other 
mixing  and  curing  additives.  With  the  exception  of  the  80nm  nano- aluminum,  all  Series  A  ingredi¬ 
ents  are  purchased  from  Firefox  Enterprises,  a  commercial  hobby  rocketry  and  fireworks  supplier. 
Ethyl  hexyl  acetate  (EH A),  polydimethylsiloxane  (PDMS),  and  lecithin  are  all  used  during  the 
mixing  process  to  allow  for  easier  mixing,  and  isonate  143L  and  HX-878  (tepanol)  are  used  for 
curing.  A  summary  of  the  ingredients  is  shown  in  Table  6. 

For  a  preliminary  test  of  the  simulation  code  over  a  range  of  compositions,  the  solids  compo¬ 
sition  was  varied  over  three  formulations.  Based  on  historical  Purdue  propellants,  the  total  solids 
composition,  coarse  to  fine  ratio  of  AP  (cAP/fAP),  and  binder  composition  were  held  constant 
at  a  nominal  76.53%,  0.648,  and  23.47%,  respectively.  Within  the  solids  content,  the  AP  percent 
composition  decreased  as  needed  while  the  nano- aluminum  loading  increased  from  nominal  values 
of  6%,  12%,  and  18%.  A  summary  of  the  actual  solids  and  binder  composition  are  listed  in  Tables 
7  and  8,  respectively. 

For  Series  A  propellants,  80nm,  uncoated  nano-aluminum  powder  from  NovaCentrix  was  chosen 
for  its  readily  available  supply  at  Purdue  University.  In  order  to  determine  its  particle  size  distri- 
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bution,  SEM  or  TEM  imaging  is  needed.  For  these  mixes  specific  imaging  is  not  available  for  the 
specific  nano-aluminum  inventory  T2R0-24-3,  but  by  argument  of  similarity  the  powder  is  assumed 
to  be  the  same  as  two  other  NovaCentrix  80nm,  uncoated  powder  inventory  lots.  Based  on  SEM 
imaging,  shown  in  Fig.  18,  and  manual  analysis,  size  distribution  of  nano-aluminum  powders  can 
be  defined  using  probability  density  and  cumulative  probability.  Plots  are  presented  in  Figs.  19(a) 
and  19(b)  for  lots  M2665B_6410  and  M2666-5406,  respectively,  and  summarized  in  Table  9. 


Figure  18:  SEM  Imaging  of  NovaCentrix  80nm  Nano-Aluminum,  Lot  M2666-5406 


Table  6:  Series  A  Propellant  Constituents 


Ingredient 

Function 

Source 

80nm  Nano- Aluminum 

Ammonium  Perchlorate  (60-130/mi  Diameter) 
Ammonium  Perchlorate  (200/mi  Diameter) 
Hydroxyl- Terminated  Polybutadiene  (HTPB) 

2  Ethyl  Hexyl  Acetate  (EHA) 

Lecithin 

Polydimethylsiloxane  (PDMS ) 

Isonate  143L  -  Modified  MDI 

HX-878  (Tepanol) 

Fuel 

Oxidizer 

Oxidizer 
Binder/Fuel 
Plasticizer 
Processing  Aid 
Anti-foam  Agent 
Curing  Agent 
Bonding  Agent/Curative 

NovaCentrix 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 
Firefox  Enterprises 

33 


Size  [nm) 


o 


n 

8 

8 

I 


(a)  Lot  M2665B-6410 


(b)  Lot  M 2666-5406 


Figure  19:  NovaCentrix  80riin  Nano- Aluminum  Size  Distribution 


Table  7:  Series  A  Solids  Composition 


Formulation 

%Nano- Aluminum 

%Coarse  AP 

%Fine  AP 

cAP/fAP  Ratio 

Total 

Al 

(5.00 

27.72 

42.81 

0.648 

76.53 

A2 

12.03 

25.43 

39.22 

0.648 

76.68 

A3 

18.00 

23.00 

35.50 

0.648 

76.50 

Table  8:  Series  A  Binder  Composition 


Formulation 

%HTPB 

%EHA 

%Lecithin 

%PDMS 

%Isonate  143L 

%Tepanol 

Total 

Al 

14.84 

4.58 

0.50 

Negligible 

2.95 

0.60 

23.46 

A2 

14.84 

4.59 

0.50 

Negligible 

2.96 

0.60 

23.32 

A3 

14.83 

4.58 

0.51 

Negligible 

2.94 

0.65 

23.50 
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Table  9:  NovaCentrix  80nm  NanoAluiniiium  Size  Distribution 


Inventory  Lot 

Mean  Size  [nm] 

Standard  Deviation  [nm] 

M2665B.6410 

77.90 

38.38 

M2666.5406 

79.14 

38.43 

Average 

78.52 

38.41 

For  larger  micrometric-size  powders  such  as  the  coarse  and  fine  ammonium  perchlorate,  size 
distribution  can  be  determined  using  a  Svmpatec  helium- neon  laser  optical  system  (HELOS),  shown 
in  Fig.  20,  and  WINDOX  4  data  analysis  software.  With  this  setup  dry  AP  powder  is  dusted  in 
front  of  a  laser  source  and  the  resulting  diffraction  pattern  is  captured  by  an  R3  or  R6  Fourier  lens 
whose  nominal  measurement  ranges  are  0.5  -  175/im  and  9-1750/mi,  respectively.  Ftom  prior  work 
with  the  R3  lens,  all  AP  powders  are  composed  of  a  very  small  amount  of  particles  less  than  10/mi 
in  diameter,  and  consequently  data  primarily  comes  from  the  R6  lens  unless  otherwise  noted.  Figs. 
21(a)  -  21(b)  and  22(a)  -  22(b)  plot  the  probability  density  and  cumulative  probability  of  Firefox 
AP  crystals;  the  mean  and  standard  deviation  are  summarized  in  Table  10. 


Figure  20:  Sympatec  HELOS  Particle  Sizer 


Table  10:  Series  A  Ammonium  Perchlorate  Size  Distribution 


Reference  Name 

Sample 

Mean  Size  [/im] 

Standard  Deviation  [/im] 

Sample  1 

167.76 

99.89 

Firefox  60-130/mi 

Sample  2 

165.03 

84.10 

Average 

166.40 

92.00 

Sample  1 

217.45 

79.41 

Firefox  200/im 

Sample  2 

206.11 

75.46 

Average 

211.78 

77.44 

With  a  baseline  propellant  matrix  established  for  steady-state  tests,  the  next  series  of  propel¬ 
lants,  Series  B  focused  on  using  higher  quality  ingredients  for  more  consistent  results  for  dynamic 
analysis.  For  the  first  composition,  propellant  Bl,  the  percent  composition  is  based  on  Galfetti 
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(b)  Sample  2 


Figure  21:  Firefox  60-130/im  Ammonium  Perchlorate  Size  Distribution 


too 

ra-  (mitronl 


**  I 


(a)  Sample  1 


(b)  Sample  2 


Figure  22:  Firefox  200/mi  Ammonium  Perchlorate  Size  Distribution 


et  al’s  baseline  formulation  [2  using  newly  acquired  NovaCentrix  SGiini  nano-aluminum  and  20/iin 
and  200 (im  AP  from  Alliant  Techsystems  Inc  (ATK).  Since  SEM  imaging  is  not  yet  available  for 
the  new  stock  of  80nm  nano- aluminum,  the  size  distribution  is  again  assumed  to  be  the  same  as  the 
data  listed  in  Table  9  since  they  are  expected  to  be  very  similar  to  other  80nm  powders.  The  AP 
size  distribution  is  available  from  the  HELOS  particle  sizer,  shown  in  Figs.  23(a)  -  23(b)  and  24(a)  - 
24(b),  where  the  R3  Fourier  lens  is  used  for  the  20/zm  powder  due  to  its  nominally  smaller  size.  The 
multiple  peaks  in  Figs.  23(a)  and  23(b)  are  due  to  moisture  exposure  during  sample  transportation 
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Probability  Density 


which  results  in  noticeable  clumping.  Consequently,  further  analysis  is  needed  before  calculating 
the  statistical  mean  and  standard  deviation  of  the  20^m  AP. 
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Figure  23:  ATK  20/ini  Ammonium  Perchlorate  Size  Distribution 
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(a)  Sample  1 


(b)  Sample  2 


Figure  24:  ATK  200/xm  Ammonium  Perchlorate  Size  Distribution 
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Cumulative  Density 


For  comparison,  the  HELOS  particle  size  data  is  compared  to  manufacturer  provided  specifi¬ 
cations  generated  either  by  sieve  based  measurements  (Rotap)  for  the  200/mi  AP  or  air  classifying 
mill  (ACM)  for  the  20/xm  AP.  Shown  in  Figs.  25(a)  and  25(b),  qualitative  analysis  suggests  that 
the  HELOS  consistently  overestimates  the  mean  particle  size  by  a  small  amount;  the  general  shape 
of  the  HELOS  curve  shows  no  correlation.  To  determine  which  measurement  method  is  most  cor¬ 
rect,  photographic  imaging  and  manual  analysis  of  particle  sizes  is  needed.  Available  statistical 
information  is  tabulated  in  Table  11  and  the  complete  propellant  B1  formulation  is  shown  in  Tables 
12  and  13. 
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(a)  ATK  20/tm  ACM  and  HELOS  (b)  ATK  200/im  Rotap  and  HELOS 


Figure  25:  Comparison  of  Manufacturer  and  Purdue  Particle  Size  Data 


Table  11:  Series  A  Ammonium  Perchlorate  Size  Distribution 


Reference  Name 

Measurement  Type 

Mean  Size  [/mi 

Standard  Deviation  [/mi] 

ATK  20/xm 

ACM 

20.11 

19.92 

HELOS 

n/a 

n/a 

ATK  200/im 

Rotap 

HELOS 

n/a 

214.29 

n/a 

87.11 

Table  12:  Series  D  Solids  Composition 


Formulation 

%Nano- Aluminum 

%Coarse  AP 

%Fine  AP 

cAP/fAP  Ratio 

Total 

B1 

15.00 

54.40 

13.60 

4.00 

83.00 

5.2  Experimental  Setup)  and  Results 

For  Series  A  propellants,  the  primary  goal  was  to  provide  a  wide  range  of  steady-state  burn  rate 
data  for  computer  model  calibration.  To  accomplish  this,  propellants  were  mixed  by  hand  and 
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Table  13:  Series  B  Binder  Composition 


Formulation 

%HTPB 

%EHA 

%Lecithin 

%PDMS 

%Isonate  143L 

%Tepanol 

Total 

Bl 

10.74 

3.31 

0.37 

Negligible 

2.14 

0.44 

17.00 

oven-cured  for  at  least  two  days.  Cast  propellant  blocks  were  then  cut  into  6x6x38mrn  (0.25  x  0.25 
x  1.5  in)  strands  and  inhibited  with  epoxy  on  all  sides  except  the  top  surface  to  provide  for  a  planar 
burn.  Although  a  strand  heater  has  been  tested  operational,  initial  tests  minimize  the  effects  of 
temperature  variation.  Room  temperature  was  kept  between  15.5°C  and  18.3°C  (60-65°F)  and 
strands  were  left  in  the  open  at  least  30  minutes  to  come  to  thermal  equilibrium.  Strands  were 
then  loaded  inside  a  windowed  pressure  vessel  and  ignited  with  nichrome  wire  at  initial  pressures  of 
34,  47.6,  68,  95.2,  or  129.3  atmospheres  (500,  700,  1000,  1400,  and  1900  psia).  A  camera  recorded 
the  burning  propellant  strand  at  30  frames  per  second  and  the  burn  surface  progression  is  manually 
tagged  frame-by-frame  using  video  analysis  tool,  Tracker  [33].  Figs.  26  and  27  demonstrate  a  typical 
planar  burn  and  show  how  Tracker  is  able  to  follow  the  burn  surface  at  each  frame  throughout  the 
burn. 


Figure  26:  Planar  Propellant  Burn 

The  data  is  then  plotted  on  a  log- log  graph  and  fitted  to  using  St.  Robert’s  Law,  r&  =  apn, 
where  r&  is  the  burn  rate  in  cm/s,  a  is  the  burn  coefficient,  p  is  the  vessel  pressure  in  atm,  and  n  is 
the  pressure  exponent.  Data  for  Series  A  propellants  is  shown  in  Fig.  28  as  well  as  a  comparison 
burn  at  68atm  with  propellant  Bl.  Series  A  burn  coefficients  and  pressure  exponents  are  listed  in 
Table  14. 


Table  14:  Series  A  Burn  Rate  Constants 


Formulation 

Burn  Coefficient,  a 

Pressure  Exponent,  n 

Curve  Fit  R2 

A1 

0.333 

0.268 

0.858 

A2 

0.300 

0.332 

0.910 

A3 

0.369 

0.382 

0.491 
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Figure  27:  Tracker  Video  Analysis  Software 
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Figure  28:  Series  A  Propellant  Burn  Rate 


Evaluating  St.  Robert’s  Law  using  the  coefficients  in  Table  14  then  allows  comparison  with  the 
theoretical  burn  rate  predicted  by  the  simulation  code  at  68atm.  The  result  is  shown  in  Table  15. 
To  explain  the  large  discrepancy  between  theoretical  and  experimental  A3  propellant,  in  comparison 
to  other  Series  A  strands,  A3  showed  significantly  more  air  pockets  which  would  consequently  result 
in  a  higher  measured  burn  rate.  Surface  scans  of  propellants  A1-A3  as  well  as  propellant  Bl  are 
presented  in  Figs.  29(a)  -  29(d)  for  comparison. 

Although  Sam  met  et.  al  [6]  have  demonstrated  little  difference  between  hand  and  mechanical 
mixes,  propellant  Bl  was  made  using  Purdue’s  Ross  vacuum  mixer  in  efforts  to  maximize  the 
consistency  and  repeatability  of  future  work.  Qualitative  analysis  of  Figs.  29(a)  -  29(d)  suggests 
slight  improvement  in  AP  crystal  dispersion,  but  markedly  minimizes  the  number  of  air  pockets. 
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Table  15:  Comparison  of  Theoretical  and  Experimental  Burn  Rate 


Formulation 

Theoretical  Burn  Rate  [em/s] 

Experimental  Burn  Rate  [em/s] 

Pereent  Difference 

A1 

1.0725 

1.032 

1.94 

A2 

1.1439 

1.218 

3.12 

A3 

1.0902 

1.849 

25.8 

(c)  A3 


(d)  B1 


Figure  29:  Comparison  of  Series  A  and  B  Propellant  Quality 


Additionally,  unlike  previous  mixes,  further  steps  were  taken  to  ensure  adequate  rheology  in  B1  by 
ultrasonieally  dispersing  the  nano-aluminum  powder  a  hexane  solution  to  reduee  agglomeration. 

To  start  the  dynamic  analysis  objective  of  the  program,  tests  for  Bl  consisted  of  cutting  pro¬ 
pellant  blocks  into  smaller  6x6x6mm  cubes  (0.25  x  0.25  x  0.25in)  and  inhibited  on  all  but  the  top 
face  with  a  thin  layer  of  polyurethane  for  better  burn  surfaee  imaging.  The  CO2  laser  was  aimed 
and  calibrated  to  output  7.5  watts  of  power  using  a  Coherent  LM-200  laser  power  detector  and  the 
ignition  transient  was  eaptured  with  a  high-speed  Phantom  eamera  recording  at  10,000  frames  per 
seeond.  Figs.  30(a)  -  30(c)  show  the  ignition  for  a  const  ant- power  laser  pulse  over  the  entire  burn, 
and  sinee  only  few  tests  were  able  to  be  completed  before  failure  of  the  refleet ive  mirror,  no  better 
images  were  taken. 

6  Future  Work 

Due  to  the  very  consistent  results  from  this  batch,  further  A-series  mixes  will  continue  to  use  the 
same  hand-mix  and  easting  procedures  to  provide  calibration  data  over  a  wide  range  of  operating 
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(a)  t  =  0 


(b)  t  =  +0. 0223s 


(c)  t  =  +0. 1223s 


Figure  30:  High  Speed  Imaging  of  Laser  Ignition 


pressures.  Once  complete,  experimental  efforts  will  move  on  to  mimicking  existing  formulations 
2  to  determine  repeatability,  while  continuing  to  perfect  casting  procedures  in  quartz  tubes  in 
preparation  for  characterization  at  high  initial  propellant  temperatures  and  dynamic  transient  ef¬ 
fects  caused  by  CO2  laser  ignition.  In  addition  to  continuing  high-speed  imaging  of  laser  ignition 
transients,  work  should  continue  with  integration  of  the  spectroscopy  system  and  plans  have  been 
made  to  develop  a  control  system  for  the  CO2  laser  to  generate  variable  pulse  widths.  Concur¬ 
rently,  progress  should  continue  on  further  defining  Series  B  propellants  and  re-evaluating  Series  A 
propellants  at  elevated  temperatures  to  determine  the  temperature  sensitivity  of  this  formulation. 
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